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Introduction 


The common feature of the network problems which we consider in this 
unit is that they can all be modelled by bipartite graphs. Some of these 
problems (such as the marriage problem introduced in Section 1) originated 
as problems in pure mathematics, whereas others (for example, the 
transportation problem discussed in Section 4) arose from practical problems 
in industry and commerce. 


Most of the problems discussed in this unit are optimization problems for 
which we give methods of solution in the form of algorithms. The three 
basic optimization problems which we discuss — the matching problem, 
the assignment problem and the transportation problem — are closely 
related to the flow problems considered in Networks 1 and can be solved by 
using the maximum flow algorithm. However, the algorithms given here 
are more efficient, since they make use of the fact that all the graphs 
involved are bipartite graphs. 


In Section 1, Matching problem, we discuss the marriage problem and some 
related problems involving matchings in bipartite graphs. 


In Section 2, Maximum matching problem, we consider problems in which 
we wish to find matchings in the corresponding bipartite graphs with as 
many edges as possible. We introduce the concept of an alternating path, 
and present an algorithm for finding a maximum matching based on this 
concept. 


Section 3, Assignment problem, is concerned with assignment problems in 
which we have costs associated with the possible assignments and we 
wish to find matchings with least overall cost. The Hungarian algorithm 
which we give for solving this problem is based on the procedure used in 
the algorithm for finding a maximum matching. We also discuss a 
variation on the basic assignment problem, known as the bottleneck 
assignment problem. 


In Section 4, Transportation problem, we discuss transportation and 
transhipment problems. The Hungarian algorithm for the transportation 
problem uses a similar procedure to the algorithm for the assignment 
problem. 


1 Matching problem 
1.1 Marriage problem 


In this section we discuss some problems involving matchings in bipartite 
graphs. We begin by presenting three examples of the type of problem 
which we consider. 


Example 1.1 


Four jobs are to be assigned to four applicants. Applicant 1 is qualified to 
do only jobs A and B; applicant 2 is qualified to do jobs A, B, C and D; 
applicant 3 is qualified to do only jobs A and B; and applicant 4 is 
qualified to do only job B. Can all four applicants be assigned to different 
jobs for which they are qualified? NO S 


Example 1.2 


Four men A;, Az, A3 and A, each know some of four women B,, B,, Bz and By. 
Man A, knows women B, and Bz; man A, knows all four women; man A; 
knows women B, and B,; and man A, knows woman By, Is it possible for all 


the men to marry women they know? \\O # 


+ 


Example 1.3 


A company has four lorries and four depots. Lorry 1 can carry loads up to 
1; tonnes; lorry 2 can carry loads up to 2 tonnes; and lorries 3 and 4 can each 
carry loads up to 1 tonne. Depot 1 needs a lorry to carry loads up to e 
tonnes; depot 2 needs a lorry to carry loads up to 1 tonne; depot 3 needs a 
lorry to carry loads up to 1; tonnes; and depot 4 needs a lorry to carry up to 2 
tonnes. Can each depot be allocated a suitable lorry? =o a 


In fact, all three examples can be expressed as exactly the same 
mathematical problem. This becomes obvious if we represent each 
situation by a graph. For Example 1.1, the following table shows the jobs 
for which each applicant is qualified. A graphical representation of this 
information is given on the right. 


applicants jobs applicants jobs 
1 A,B 1 A 
2 A,B,C, D Z B 
3 A,B a a 
4 B 4 D 


In the above graph, the black vertices represent applicants and the white 
vertices represent jobs. A black vertex is joined by an edge to a white 
vertex if the corresponding applicant is qualified to do the corresponding 
job. Since only vertices of different colours are joined, the graph is a 
bipartite graph. 


If we represent Examples 1.2 and 1.3 by graphs in a similar way, we find 
that the graphs are identical to the one obtained above. For Example 1.2, 
we simply substitute ‘men’ for ‘applicants’, and ‘women’ for ‘jobs’. For 
Example 1.3, we substitute ‘depots’ for ‘applicants’, and ‘lorries’ for ‘jobs’. 


Since the three problems are essentially the same, we restrict our 
discussion to the solution of Example 1.1. 


Assigning an applicant to a job corresponds to choosing an edge of the 
bipartite graph. We can easily assign three of the applicants to different 
jobs; for example, we can assign applicant 1 to job A, applicant 2 to job D, 
and applicant 3 to job B. 


This arrangement is shown in the following graph, where thick edges 
represent assignments of applicants to jobs. 


applicants jobs 


1 A 
2 B 
3 c 
4 D 


The pairing of some or all of the elements of one set with elements of a 
second set is cailed a matching. The matching represented by the thick 
lines of the above graph comprises the following pairings: 


applicant 1 with job A 
applicant 2 with job D 
applicant 3 with job B 


A matching in a bipartite graph corresponds to a set of edges no two of 
which meet at the same vertex. The matching represented by the above 
graph contains three such edges, as shown by the thick lines. 


W\ 


As stated, each problem is an existence 
problem; the actual assignment is a 
construction problem. 


A formal definition of a matching ina 
graph is given in Graphs 3. In this unit 
we are interested in matchings in 
bipartite graphs in which each edge of 
the matching joins a vertex in one set 
to a vertex in the other. 


We now see that our original question: 


can all four applicants be assigned to different jobs for which they 
are qualified? 


is equivalent to the question: 

can we find a matching containing four edges in the bipartite graph? 
In fact, it is not possible to find such a matching in this graph. We say 
that the matching with three edges is a maximum matching — that is, it 
is a matching with the largest possible number of edges. 
Problem 1.1 


Suppose that four applicants apply for four jobs and that each applicant 
is qualified for some of these jobs, as shown in the following table. 


applicants jobs 


: A, 
Z A,B 
3 A; 
4 B 


(a) Draw the bipartite graph which represents this situation. 


(b) By finding a maximum matching in the bipartite graph, find the 
maximum number of applicants that can be assigned to different jobs 
for which they are qualified. 


In the above problem, it is possible for all four applicants to be assigned to 
jobs for which they are qualified. This is in contrast to our previous 
example, where at most three of the four workers can be assigned to 
different tasks for which they are qualified. For the general problem in 
which we are trying to allocate a set of applicants to various jobs, we 
should like to know under what conditions it is possible for all the 
applicants to be assigned to jobs for which they are qualified. 


This problem was formulated in 1935 as a problem in group theory. It was 
later restated in terms of marriages, as in Example 1.2. For this reason, it 
is often called the marriage problem. We shall follow this tradition of 
discussing matching problems in bipartite graphs in terms of marriages, 
rather than in terms of more practical problems such as the assignment of 
workers to tasks or applicants to jobs. This has the advantage that it 
helps us to keep clear the distinction between matching problems and the 
related, but more general, assignment problems which we discuss in 
Section 3. 


We state the marriage problem formally as follows. 


Marriage problem 


Given a set of men, each of whom knows some women from a given set of 
women, under what conditions is it possible for all the men to marry 
women they know? 


One condition which must obviously be satisfied is that the total number 
of women must be equal to or greater than the total number of men. We can 
go further than this, and state that any two men must collectively know at 
least two women. More generally, for each subset of m men, the m men 
must collectively know at least m women, for all possible values of m. 


It is assumed that each man can 
marry only one woman, and each 
woman can marry only one man. 


This is a necessary condition for it to be possible for all the men to marry 
women they know. What is not so obvious is that it is also a sufficient 
condition. This result is known as the marriage theorem. 


Theorem 1.1: marriage theorem This theorem was proved (in a 


different formulation) by Philip Hall 


A necessary and sufficient condition for there to be a solution to the | ‘ —. 
In . 


marriage problem (that is, for every man to be able to marry a woman 
he knows) is: 


for every subset of m men, the m men collectively know at least 
m women, for all values of m in the range 1 < m <n, where 7 is the 
total number of men. 


We prove this theorem at the end of this section. First, we examine how it 
applies to some specific cases. 


We begin by applying the marriage theorem to Example 1.2. The 
following table shows which women are known to each of the four men. 


men women known 
Ay B,, By 

Az B,, By, B3, Bg 
Ag B,, By 

Ay Bo 


The marriage theorem states that for all the men to be able to marry 
women they know, the men in every subset of {A;, Az, A3, Ag} of size m 
must collectively know at least m women, for each possible value of m. 
This condition is known as the marriage condition. 


So, to apply the condition systematically, we first list all possible subsets 
of men. 


There are four men, so the total number of subsets is 2*. However, we need The number of subsets of a set of 
look at only 2* — 1 = 15 subsets, since we obviously do not need to consider 1 elements is 2”. 
the empty set. 


The four men collectively know four women, so the marriage condition is 
satisfied for the subset {A,, A2, A3, Ag}. The following table shows 
whether the marriage condition is satisfied for each subset of three men. 


numberof Is marriage 


number of women condition 
subset of men women known men m known satisfied? 
{A,, Az, A3} B,, Bo, Bs, Bg 2 . yes 
{A,, Az, Ag} B,, Bo, Bz, Bg 2S 4 yes 
{A,,A3, Ag} B,, Bo 3 2 no 
{Az, Az, Aq} B,, Bo, Bs, Bg 3 + yes 


For the third of these subsets, the three men concerned collectively know 
only two women, so the marriage condition is not satisfied, and therefore 
all four men cannot marry women they know. For this example, we do not 
need to examine any more subsets, but if the marriage condition had been 
satisfied for all four subsets of three men, then we should have had to 
continue with the subsets containing two men, and if necessary, the subsets 
containing only one man. 


Problem 1.2 


Suppose that each of the three men A,;, A7 and A3 know some of the four 
women B,, By, Bz; and By, as shown in the following table. 


men women known 
Ay B,, Bz, Bg 

Ay By 

Ag Bo, Bs 


(a) Construct a table (similar to the one for the example in the text) 
showing for each possible subset of men, the number of men in that 
subset and the number of women the men collectively know. For each 
subset, indicate whether the marriage condition is satisfied. Use the 
table to decide whether it is possible for all three men to marry 
women they know. 


(b) If Gin the general case) there are n men, how many subsets of men do 
you need to consider when applying the marriage theorem? Check 
that your answer for the case n = 3 agrees with the number of entries 
in your table. 


(c) Draw the bipartite graph representing this problem. Find a 
maximum matching in this graph corresponding to the maximum 
number of men who can marry women they know. 


Problem 1.3 


Consider four men A;, Az, Az and A,and four women B,, B,, B3and By, The 
following table shows which women are known to each of the four men. 


men women known 
Ay By 

A» By, Bs 

A3 B,, By, By 

Ag B3 


(a) Draw the bipartite graph representing this situation. 


(b) Write down the subsets of men (if any) for which the marriage 
condition is not satisfied. 


(c) By inspection, find a maximum matching in the bipartite graph, and 
hence find the maximum number of men who can marry women they 
know. " | 


So far we have been mainly concerned with the question posed by the 
marriage problem: 


under what condition is it possible for all the men to marry women 
they know? 


We have seen that this question is answered by the marriage theorem. 
When it is not possible for all the men to marry, it is natural to ask a 
second question, called the modified marriage problem: 
what is the maximum number of men who can marry women they 
know? 
In a practical version of this type of problem (such as assigning workers to 
tasks), we should probably be most interested in actually constructing a 
maximum matching. We express this problem as a third question: 


how can we determine which man should marry which woman so 


that the maximum possible number of marriages can take place? Our answer to the third question 
In th b hall b = A aad: d takes the form of an algorithm. 
n the next subsection we shall be mainly concerned with the secon We describe such an algorithm 


question. in Section 2. 


1.2 Modified marriage problem “Ase. 


In the solution to Problem 1.3, we found by trial and error that of the four 
men A;, Az, A3and A,, only three can marry women they know. We now 
show how to prove that the maximum number of men who can marry is 
three, without having to construct a possible scheme of marriages. 


For convenience, we set out below a table and a bipartite graph showing 
which men know which women. 


men women known men women 


A; By Ay By 
A By, Bg A By 
A3 By, Bo, By Ag B; 
Ay B, Ag B, 


We found in Problem 1.3 that the marriage condition is satisfied for all 
subsets of men except for the subset containing the three men A;, A and Ay, 
who collectively know only two women, B, and B;. By the marriage 
theorem, it is not possible for all four of the men to marry. 


Since the marriage condition breaks down because the three men Aj, A, 
and A, know only two women, it seems likely that if we remove one of 
these men, then the marriage condition will be satisfied for all possible 
subsets of the remaining three men. This is indeed the case. For example, 
if we remove Aj, then the subset {A,, A, A,} becomes the subset {A,, A,} 
and A, and A, collectively know two women. 


men women 

B, 

a B, 
A3 Ba 

Ay By 


All the other subsets which contained A, and which previously satisfied 
the marriage condition still satisfy the marriage condition, so all three 
men can marry women they know. Although this is fairly obvious if we 
remove A,, it is certainly not obvious that we can apply similar reasoning 
to a more complicated example. | 


Let us try another approach. Faced with a situation where the marriage 
condition was not satisfied, we converted the situation to one where the 
- Marriage condition was satisfied, by removing one of the men. However; 
we can also satisfy the marriage condition for all subsets of men by adding 
an extra women B;, who is known by all the men. The men in the subset 
{A,, Az, Ag} now know three women B,, Bz and Bs, and all the other subsets 
of men now know one more woman than before. Therefore, by the marriage 
theorem, all four men can marry women they know. If we now remove the 
extra woman Bs, we are left with three men who are married and one (the 
one who was married to Bs) unmarried. We have thus shown that it is 
possible for three of the men to marry women they know. 


We now turn our attention to the general problem of finding the maximum 
number of men who can marry when the marriage condition is not 
satisfied. 


If there are n men, then the marriage condition states that they can all 
marry women they know if 


(number of women known by each subset of m men) 2 m, 
for all values of m in the range 1 <m <n. 


If the marriage condition is not satisfied, we can write instead that 
n —d men can marry women they know if 


(number of women known by each subset of m men) 2m — d, 
for all values of m in the range 1 <m <n, and for some positive 
integer d. 


This is called the modified marriage condition. We prove this result by an 
argument similar to the one we used above for the example of Problem 1.3. 
All we have to do to prove the result is to add a group of d women who are 
known by all the men. With the enlarged group of women, the marriage 
condition is satisfied, so all the men can marry. Removal of the extra d 
women leaves n —d men who are married, and d men who are not married. In 
other words, in the original problem, n — d men can marry. 


If we put r = n-—d, we can write the modified marriage condition as: 


(number of women known by each subset of m men) 2m-—(n—-71), 
for all values of m in the range 1 <m <n. 


The maximum value of r for which this modified marriage condition is 
satisfied is the largest number of men who can marry women they know. We 
can rearrange the inequality in this marriage condition so that only r 
appears on the right-hand side: 


(the number of women known by each subset of m men) + (1 —m) 21, 
for all values of m in the range 1 <m <n. 


It is clear that the maximum value of r for which this condition is satisfied 
is equal to the minimum value of the expression on the left-hand side over 
all subsets of m men, over all values of m. We summarize this result in the 
following theorem. 


Theorem 1.2: modified marriage theorem 


If a group of n men each know some of a group of women, the maximum 
number of men who can marry women they know is equal to the minimum 
value of the expression: 


(number of women p known by a subset of m men) + (n - m), 
for any subset of m men, and for all values of m in the range 
LSosn. 


Worked problem 

Suppose that five men A, Az, A3, Agand A; know four women By, By, Bz and 
B,, as shown in the following table and bipartite graph. What is the 
maximum number of men who can marry women they know? 


men women known men women 
Ay 8.8; A; By 
A By, B3 Ag By 
A3 By, By A3 B; 
A, By, By, Bs, By Ag By 
As B, | As 


Solution 


To answer this question, we need to find the minimum value of the 
expression in the modified marriage theorem. To do this, we need consider 
only those subsets of men for which the original marriage condition is not 
satisfied. By inspection of the table or the graph, we see that there are 
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Note that if we put d = 0, then r is 
equal to n, and this condition reduces 
to the original marriage condition. 


This theorem is a minimax theorem, like 
the max-flow min-cut theorem and 
Menger's theorems in Networks 1. 


In stating this theorem, we assume that 
the original marriage condition is not 
satisfied. Otherwise the minimum value 
of the expression may be greater than n. 


only two such subsets: {A;, Az, A3, As} and {A;, Az A3, Ag, As}. We draw up 
the following table. 


number of 
numberof © women 
subset of men women known menm knownp p+(n-m) 
{A,, Az, A3, As} B,, By, B3 4 s 4 
{A,,Az,A3,A4,As5} By, By, B3, Bg 5 d 4 


For both subsets, the value of the expression p + (n - m) is 4. Hence exactly 
four men can marry women they know. = 


In the above worked problem, we found that at most four men can marry. 
The next question which arises is: which four of the five men can marry 
women they know? Or, to put it another way, which one man should we 
remove so that the remaining four can marry? To answer this question, we 
look at the two subsets for which the marriage condition is not satisfied. 
We must remove a man who is in both of these subsets. Clearly, there are 
four possibilities: A;, Az, Az or As. We remove each of A;, Az, Aj and As in 
turn, and find a marriage scheme for the remaining four men in each case. 
All such marriage schemes involving the maximum number of men (four) — 
are shown below. 


B, Ay B, Ay B, A, 
A B, By A By Ay 
A3 B; A3 Bz Bz A3 
A, B, A, By Ay By A, 
As As As 
A, removed A, removed Az removed 


If A, is removed, then it is obvious that the remaining four men cannot all 
marry women they know, since none of them knows By. 


Problem 1.4 


Suppose that five men A;, Az, A3, Ag and A; know four women B;, Bz, B3 
and By, as follows. 


men women known 
Ay By, Bs 

Az Bo, Bz 

A By 

Ag B3 

As By, Bo, Bg 


(a) Show that the marriage condition is not satisfied. 
(b) Draw a bipartite graph to represent this problem. 


(c) List all the subsets of men for which the marriage condition is not 
satisfied. 
For each of these subsets, find the value of the expression 
p + (n—™m) in the modified marriage theorem. 


(d) What is the largest number of men who can marry women they know? 
List all the subsets of men containing this number of men who can all 
_ marry women they know. 


By Ay 
B, Ay 
B; Ag 
By Ag 


As removed (two possibilities) 


One way of doing this is to construct a 
table as in the solution to Problem 1.2. 


Alternatively, you may find the 


required subsets by inspection, using 


the graph representation. 
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We have now answered the first and second of the three questions we 
listed earlier. The remaining question is: 


how can we determine which man should marry which woman so 
that the maximum possible number of marriages can take place? 


To answer this question, we might try to adapt the procedure based on the 
modified marriage theorem which we used in the solution to the previous 
problem. This procedure can be converted into an algorithm fairly easily, 
but such an algorithm is not very suitable for solving a practical form of 
the marriage problem in which, for example, we wish to assign a number 
of workers to tasks for which they are qualified. The method is unsuitable 
for the following reasons. 


e Although the procedure tells us which subsets of men can all marry 
(or which workers can be assigned to tasks), it does not tell us which 
man should marry which woman. For small groups of men it is a 
fairly simple matter to construct a suitable marriage scheme, but for 
larger numbers we need a systematic method. 


e — The procedure is an exhaustion method which involves looking at 
all possible non-empty subsets of the set of nm men. The number of such 
subsets is 2” — 1, which increases very rapidly as n increases. The 
time taken to perform an algorithm based on this procedure is 
therefore roughly proportional to 2”, that is, to an exponential 
function of n. Since the time taken increases so rapidly with 1, it 
would not be practicable to use the procedure for problems involving 
large values of n, even when using a computer. 


e The procedure is wasteful, since it effectively carries out a search of 
all possibilities and finds all the subsets of men which contain the 
largest number who can marry (or the largest number of workers who 
can all be assigned to different tasks). In a practical situation, we 
normally require only one solution. 


What we need is an efficient algorithm which constructs just one maximum 
matching. We describe such an algorithm in Section 2. Before doing this, 
we give two different proofs of the marriage theorem. 


Not 


1.3 Two proofs of the marriage theorem j.2~<4 


For convenience, we repeat the statement of the marriage theorem. 


Theorem 1.1: marriage theorem 


A necessary and sufficient condition for there to be a solution to the 
marriage problem (that is, for every man to be able to marry a woman 
he knows) is: 


for every subset of m men, the m men collectively know at least 
m women, for all values of m in the range 1 < m <n, where n is the 
total number of men. 


Since it is obvious that the condition given in the statement of the 
marriage problem is necessary, we need prove only that it is sufficient. We 
prove this first by mathematical induction on n and then by a constructive 
proof. 


Proof by mathematical induction 
We carry out two steps. 


STEP 1 Verify that the theorem is true when n = 1. 


In this case, there is one man who knows at least one woman, so obviously 
the marriage problem has a solution. 
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By a constructive proof we mean a 
proof which actually shows how to 
construct a solution to a problem, as 
contrasted with a proof (such as a 
proof by induction) which establishes 
that a solution exists, but does not 
provide a procedure for finding one. 


STEP 2 Show that, if the theorem is true for every number of men less 
than n, then it must also be true for n men. 


We assume that the theorem is true for every number of men less than n. We 
consider n men. There are two cases which together cover all possibilities. 


Case (a) Suppose that for all m < nevery set of m men collectively know 
at least m + 1 women. 


In this case, the marriage condition is satisfied with ‘at least one woman to 
spare’ for every set of m men. We can therefore take any man and marry 
him to any woman he knows, and the marriage condition will still be 
satisfied for the remaining n — 1 men. By our assumption, the theorem is true 
for every number of men less than n, so it must be true for these 
n-imen. We can therefore marry off the n - 1 men appropriately. We 
have now married off all n men, so this completes Step 2 in this case. 


Case (b) Suppose that m <n and that there is at least one set of m men who © 


collectively know exactly m women. 


From our assumption that the theorem is true for every number of men less 
than n, it follows immediately that we can marry these m men to the 
m women, leaving n —- m men and (since we are also assuming that the 
marriage condition is satisfied) at least n -m women. 


Now any collection of r men from these n —-m men must collectively know at 
least r women; otherwise these r men, together with the set of m men, would 
collectively know less than r + m women, which is contrary to the marriage 
condition. It follows that, since we are assuming that the marriage theorem 
is true for fewer than n men, we can marry the remaining n —m men to women 
they know. Thus, if the statement is true for all numbers of men less than 1, 
then it is true for m men. This completes Step 2 in this case. 


Therefore, by the principle of mathematical induction, the theorem is 
true for all positive integers n. a 


Constructive proof 


To prove that the marriage condition is a sufficient condition for there to be 
a solution to the marriage problem, we show that, starting with m men 
paired with m women (where m < n), it is possible to increase this to a 
pairing of m+ 1 men. There is no difficulty in finding an initial pairing, 
since we can if necessary start with m = 1; that is, by pairing any man with 
a woman he knows. We can then successively increase the number of men 
paired until all n men are paired with women they know. 


We begin with m men paired with m women whom they know. 


If there is a man left who knows a woman not already paired, then an 
(m + 1)th pairing is immediately possible. 


Otherwise, we proceed as follows. 


We choose any unpaired man Ag. Since the marriage condition is satisfied, 
there must be at least one woman whom he knows — we choose one of these 
women B,. Suppose that B, is paired with a man Aj,. This pairing is 
indicated on the diagram by a thick line, and the fact that Ag knows B, is 
indicated by the thin line joining the corresponding vertices. 


By the marriage condition, Ap and A; must collectively know at least two 
women, namely, B, and at least one other woman B,. If Bz is not already 
paired, then A, must know By), so we pair Ag with B, and A, with B, 
(removing the original pairing of A; with B,); we thus obtain one more 
pairing, as required. 


m women 


o © 


m — 1men<x< e 
$ re) 


mmen m+1women 


eee meng 

m men m women 
a 
we O 


n-m men 


at least 
s o ( n—-m women 

fe) 

men women 
: ae 

Ay B, 
eee... 
Ay Sree + 

OBr 44 
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If By is paired with Az, then the men Ay, A; and A, must collectively know 
a third woman B3. If B3; is paired with A3, we continue until an unpaired 
woman B,,; is reached. A possible result of this procedure is illustrated on 
the left below. 


men women men women 
Ag Ag 
Ay B 1 Ay B 1 


r pairings r+1 pairings 


Note that each woman B; involved in this process must be known by at 
least one man Aj, where j < i. This fact enables us to rearrange the pairings 
as follows. First we pair woman B,,, with a man A; who knows her (where 
i<r + 1). This frees woman B;, who can be paired with some man A; 
(where j < i). We continue in this way until a woman is free who can be 
paired with Ag. This process corresponds to tracing back a path which 
starts from B,,, and proceeds by alternate thin and thick lines back to Ag. 


We have now replaced the r original pairings in the above graph with 
r+ 1 pairings. Together with the pairings not involved in this process, we 
now have a total of m + 1 pairings. : 


If there are still any men who are not paired (that is, if m + 1<n), we 
continue with this process until all m men are paired with women they 
can marry. : ra @ 


The constructive process of the second proof provides a method for finding 
a matching of the n men to n of the women. The algorithm which we 
present for finding a maximum matching is based on the procedure of this 
proof. | 


2 Maximum matching problem 


We now tackle the problem: 


how can we find a maximum matching in a bipartite graph? 


In this section we give an algorithm for finding such a matching. This 
algorithm is based on the concept of an alternating path. 
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2.1 Alternating paths 


In the following bipartite graph, a matching M consists of the two edges 
XY. and x3y3 (indicated by thick lines). The diagram on the right 
illustrates an example of an alternating path with respect to M — only 
alternate edges of this path x4 y3 x3 Y2%1 y are edges of the matching M. 


xX Y x Y There are three other alternating 
xy Y} xy Yt paths — the paths x9 y2%1¥),%2y1 and 
X4Y4- 
x2 Y2 X72 @ Y2 
x3 Y¥3 X3 ¥3 
X4 Y4 X4 OY4 
bipartite graph an alternating path 


Definition 


Let G be a bipartite graph in which the set of vertices is divided into 
two disjoint subsets X and Y. An alternating path with respect to a 
matching M in G is a path which satisfies the conditions: 


(a) the path joins a vertex x in X to a vertex y in Y; 

(b) the initial and final vertices x and y are not incident with an edge 
in M; 

(c) alternate edges of the path are in M, and the other edges are not 
in M. 


It follows that the Ist, 3rd, 5th, ... edges are not in M, and the 2nd, 4th, 
6th, ... edges are in M. Note that an alternating path has odd length. 


Problem 2.1 
Consider the following bipartite graph and matching M. 


x4 y ie 


Which of the following paths are alternating paths with respect to M? 
(a) X3Y2%y = (BD) YaXay2%YWi%1 == (C)_ Xa Ya X4Ys 


Problem 2.2 


Find all the alternating paths with respect to the matching shown in the 
following bipartite graph. 


x4 Yj 
X2 Y2 
- +s Y¥3 
X4 Ya 
X5 Y5 
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To see why the concept of an alternating path is useful for finding a 
maximum matching, consider the following diagrams. 


a Y1 *4 yi 

“7 Y2 x2 @ Y2 

se. ¥3 a“ a 

X4 OY4 X4 OY 
alternating path new matching M’ 

X4 W3 XZ Yo X1 Vy with 3 edges: 
with respect to matching X1 Yi, X3 Yo, Xa Y3 

M (with 2 edges) 


The diagram on the left shows our earlier example of an alternating path 
with respect to a matching M. In the diagram on the right, we have 
replaced the thick edges of the path by thin edges, and vice versa. In this 
way we have constructed a new matching M’ consisting of three edges — 
one more edge than in the original matching M. It follows from the 
definition of an alternating path that if we can find an alternating path 
with respect to a matching M, then we can always construct a matching M’ 
which has one more edge than M. The properties of an alternating path 
which ensure that this is the case are: 


e for an alternating path with respect to a matching M, the number of 
edges of the path not in M is always one more than the number of 
edges of the path in M; 


e the initial and final vertices of an alternating path are not incident 
with an edge in M. 


In our example, the alternating path P includes all the edges of M. If the 
alternating path P does not include all the edges of M, the new matching 
M’ consists of the edges of the alternating path P not in the matching M 
and the edges of the matching M not in the alternating path P. 


For example, if we consider the alternating path x2 yx, y, with respect to 
the matching M in our bipartite graph, we obtain an improved matching 


as follows. 
x9 Y2 xX Y2 


eae es ey aS 
X4@ OV, X4@ OV 
alternating path new matching M’ 
Xo Yo X11 with 3 edges: 
with respect to matching X1Y1, X2 Yo X3 Y3 
M (with 2 edges) 


A formal statement of this procedure is given below. 


To find a new matching with one more edge 
Find an alternating path P with respect to an existing matching M. 


Form a new matching M’ with the following edges: 


e the edges of the alternating path P not in the matching M 


and 


e the edges of the matching M not in the alternating path P. 
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2.2 Maximum matching algorithm 


This method of using an alternating path to form a new matching is the 
basis of the following algorithm for finding a maximum matching. 


To find a maximum matching, we start with a matching M, look for an 
alternating path P, and, if one is found, use it to construct a new matching 
with one more edge than before. We repeat the procedure until no more 
alternating paths can be found; the current matching is then a maximum 
matching. 


We give a formal statement of the maximum matching algorithm and 
then work through an example in the first audio-tape session. 


The algorithm is in two parts. 
Part A is the labelling procedure for finding an alternating path. 


Part B is the matching improvement procedure for finding a new matching 
with one more edge. 


In the following audio-tape session, we illustrate the use of this formal 
algorithm by working through an example. Before listening to the tape, 
read through the statement of the algorithm given here. Have this 
statement and the Audio-tape Notes to hand when you listen to the tape. 


This algorithm, in common with others presented in this unit, is written in 
a form suitable for implementing as a computer program. Its advantages 
may not be apparent from the small examples given in the tape, where 
the solutions can be quickly found by trial and error, but its usefulness is 
that it is a precise and systematic method which can be used to solve large 
examples with the aid of a computer. 


Algorithm for finding a maximum matching 


START with any matching M in a bipartite graph whose sets of vertices There is no difficulty in finding an 


aoe = (x1, %;...X,} and Y ={ypts,.-.:2Ym!- initial matching M. If necessary, it can 
consist of just one edge. 


Part A: labelling procedure 
STEP 1 Label with (+) each vertex in X which is not incident with any 
edge in the current matching M. 


If no vertex in X can be labelled, STOP: the current matching is a 
maximum matching. 


Otherwise, go to Step 2. 


STEP 2 Choose a newly labelled vertex in X, say x;, and label with (x;) 
all the unlabelled vertices in Y joined to x; by an edge NOT IN M. 


Repeat this for all newly labelled vertices in X, and then go to 
otep 3. 


STEP 3 Choose a newly labelled vertex in Y, say y;, and label with (yj) 
all unlabelled vertices in X joined to y; by an edge IN M. 


Repeat this for all newly labelled vertices in Y. 


Repeat Steps 2 and 3 until 


EITHER a vertex in Y which is not incident with an edge in M is labelled 
(this is called breakthrough), in which case go to Part B, Step 4, 


OR no such vertex exists, in which case STOP: the current matching 
is a maximum matching. 
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Part B: matching improvement procedure 


STEP 4 Find an alternating path as follows. 


Start at the breakthrough vertex, and go to the vertex indicated 
by its label. From this vertex, go to the vertex indicated by its 
label, and so on, until a vertex labelled (*) is reached. 


This path P is an alternating path. 


STEP 5 Forma new matching from: 


e the edges of the current matching M which are NOT IN the 
alternating path P; = 


e the edges of the alternating path P which are NOT IN the 
current matching M. 


Remove all labels and return to Part A, Step 1 to see whether the 
new matching can be improved further. 


€9 fa & 


2.3 Maximum matchings and maximum flows 


The problem of finding a maximum matching in a bipartite graph is 
closely related to the problem of finding a maximum flow in a network. In 
fact, a problem of finding a maximum matching can always be formulated 
as a maximum flow problem. To see how this is done, consider the 
following diagram. 


1 ~ i 
aA ~*~ 
KATE EES 
fn ee ag oa 
ay et 
a 


/ 


= 
mn 
\ 


\ 


The bipartite graph on the left is converted to a basic network as follows: 
a add a source S and a sink T; 

e replace each edge x;y; by an arc of capacity 1 directed from x; to yj; 

» for each of the vertices x; of X, add an arc Sx; with capacity 1; 

e for each of the vertices y; of Y, add an arc y;T with capacity 1. 


Suppose that we use the maximum flow algorithm to find a maximum flow 
in this network, as described in Networks 1. Then the flow along each arc 
will have the value 0 or 1. 


Furthermore, no two of the arcs shown as solid lines which then have a 
flow of 1 will meet at the same vertex, since the flow into any vertex of X 
or out of any vertex of Y cannot be greater than 1. It follows that, when the 
procedure of the algorithm has been carried out, the solid arcs with a flow 
of 1 in the network correspond to the edges of a maximum matching in the 
original bipartite graph. 


In other words, we can solve a maximum matching problem on a bipartite 
graph by converting the bipartite graph to a basic network, and using the 
maximum flow algorithm to find a maximum flow. In fact, the algorithm 
we gave above for finding a maximum matching is a special case of the 
maximum flow algorithm. 
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In the diagram, the numbers on the 
arcs are capacities, not flows. 


Problem 2.3 


In the above diagram, find (by inspection) a maximum matching in the 
bipartite graph and a corresponding maximum flow in the network. 
Which feature of the graph corresponds to the capacity of a minimum cut 
in the network? 


We can also obtain the marriage theorem by a similar construction and 
using the max-flow min-cut theorem of Networks 1. 


2.4 Computer activities 


The computer activities for this section are given in the Computer 
Activities Booklet. 


3 Assignment problem 


In the previous section we were concerned with finding a matching with 
the maximum possible number of edges in a bipartite graph. Often, there 
are several such maximum matchings. Sometimes we are concerned with 
finding just one maximum matching, and we have no reason to prefer one 
maximum matching to any other. However, if the bipartite graph 
represents a practical problem, some of the matchings may represent more 
desirable solutions to the problem than others. We can express this fact by 
assigning costs to the edges of the bipartite graph — we then look for a 
maximum matching with the lowest total cost. This type of problem is 
called an assignment problem. Let us look at an example. 


Example 3.1 


- Three people apply for three jobs. The ability of each applicant for each 
job is represented by a cost — the lower the cost, the better the ability of 
the candidate to do the job. We can represent this situation by the 
following bipartite graph. 


jobs 
Yi Y2 Y3 
x4 3 4 Z 


applicants x, a =. 


xX3 2 1 5 


cost matrix 


bipartite graph 


We have written the cost associated with the assignment of each 
applicant to each job next to the corresponding edge. Thus, for example, 


In Section 4 we shall have weights on 


both edges and vertices. To 


distinguish between the two, we refer 


to the weights on edges as costs. 
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the cost associated with assigning the applicant x; to the job y2 is 4. For 
larger graphs it is more convenient to display these costs in the form of a 
cost matrix. The cost matrix for this example is shown to the right of the 
bipartite graph. = 


The problem is to find a maximum matching for which the sum of the costs 
on the edges of the matching is a minimum. One way of proceeding would 
be to use the exhaustion method — to find all possible maximum 
matchings and then to choose one with minimum total cost. However, this 
is not an efficient procedure: for a graph with a large number of vertices, it 
would take a long time, even using a computer. What is needed is an 
efficient algorithm. Such an algorithm — called the Hungarian algorithm 
— is given below. 


In the above example, the number of people is equal to the number of jobs, 
and the graph is a complete bipartite graph. Other assignment problems 
may involve incomplete graphs or different numbers of vertices in the two 
sets. The algorithm we give can accommodate such problems, although it 
may not be possible to assign all the people to jobs. 


The algorithm is called ‘Hungarian’ in 
recognition of the Hungarian 
mathematician E. Egervary who 
introduced the idea in 1931. 


3.1 Hungarian algorithm for the assignment problem 


The algorithm starts with no initial matching in the bipartite graph 


containing the sets of vertices X = {x1,%2,...,Xn} and Y = {yj, yo, ..., Yn}. It Note that for this algorithm, the sets 


selects a set of edges of the graph in such a way that a matching involving 
these edges is at minimum cost; that is, it is not possible to find a matching 
of the same size involving other edges of the graph with a lower total 
cost. The graph containing the original vertices, together with this set of 
edges, is called a partial graph. If the matching involving the largest 
possible number of edges in the partial graph is not a maximum matching 
of the original graph, then further edges are selected and a larger 
matching is found. This process is repeated until a maximum matching at 
minimum cost is found. 


The algorithm is in three parts. 
First, the first revised cost matrix and first partial graph are constructed. 


Parts A and B are essentially the same as in the algorithm for finding a 
maximum matching given in Section 2. 


In the following audio-tape session, we illustrate the use of this 
algorithm by working through an example. Before listening to the tape, 
read through the formal statement of the algorithm given here. Have 
this statement and the Audio-tape Notes to hand when you listen to the 
tape. 


Hungarian algorithm for the assignment problem 


START with no initial matching in a bipartite graph whose sets of 
vertices are X = {X1,X2,...,X,} and Y = {y1, 2, ..., Yn}. 


Construction of first partial graph 


STEP 0a Assign weights to the vertices as follows. 


To each vertex of X, assign a weight equal to the lowest cost on 
any edge incident to that vertex. 


Decrease the cost on each edge by the weight on the vertex to 
which it is incident. 


Repeat this procedure for each vertex in Y. 


STEP 0b Construct a partial graph consisting of the vertices of the 
original bipartite graph, together with only those edges whose 
current cost is zero. 
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X and Y have the same number of 
elements. 


You can think of an original edge cost 
as being the total cost of a trip from 
one vertex to the other. We are 
dividing this cost into three parts: the 
cost of leaving the vertex in X; the cost 
of the journey (the new edge cost); 
and the cost of arriving at the vertex 
in Y. 

When Step 0a has been carried out, 
for any edge, the sum of the weights 
on the end vertices plus the new cost 
is equal to the original cost. 


Part A: labelling procedure 


STEP 1 Label with (*) each vertex in X which is not incident with any Note that the first time we meet 
edge in the current matching M. ia 1 the current matching has no 

es. 

If no such vertex exists, STOP: the current assignment is an 5 

optimum assignment. 

Otherwise, go to Step 2. 


STEP 2 Choose a newly labelled vertex in X, say x;, and label with (x;) all 
the unlabelled vertices in Y joined to x; by an edge NOT IN M. 


Repeat this for all newly labelled vertices in X. 

If no vertex in Y is labelled, STOP: the current assignment is an 
optimum assignment. 

Otherwise, go to Step 3. 


STEP 3 Choose a newly labelled vertex in Y, say y;, and label with (y;) all 
unlabelled vertices in X joined to y; by an edge IN M. 


Repeat this for all newly labelled vertices in Y. 
Repeat Steps 2 and 3 until 
EITHER a vertex in Y which is not incident with an edge in M is labelled 
(this is called breakthrough), in which case go to Part B, Step 4, 
OR no more labelling is possible, in which case go to Part C, Step 6. 


Part B: matching improvement procedure 
STEP 4 Find an alternating path as follows. 


Start at the breakthrough vertex, and go to the vertex indicated 
by its label. From this vertex, go to the vertex indicated by its 
label, and so on, until a vertex labelled (*) is reached. 


This path P is an alternating path. 
STEP 5 Formanew matching from: 


e the edges of the current matching M which are NOT IN 
the alternating path P; 


e the edges of the alternating path P which are NOT IN the 
current matching M. 


Remove all labels and return to Part A, Step 1 to see whether the 
new matching can be improved further. 


Part C: modification of partial graph procedure 


STEP 6 Find the lowest current cost on any edge of the original bipartite 

graph which 

e starts ata LABELLED vertex in X; 

e ends at an UNLABELLED vertex in Y. 

Call this lowest cost 6. 
STEP 7 (a) Increase the weight on each labelled vertex in X by 6. Note that 6 is always positive. Also, 
the total cost of a ‘trip’ from a vertex in 
X to a vertex in Y via any edge, that is, 
(c) For each edge of the original bipartite graph which joins a the sum of the two vertex weights and 


labelled vertex in X to an unlabelled vertex in Y, decrease the edge cost, is kept constant 
the current cost by 6. throughout Step 7. 


(d) For each edge of the original bipartite graph which joins an 
unlabelled vertex in X to a labelled vertex in Y, increase the 
current cost by 6. 


(b) Decrease the weight on each labelled vertex in Y by 6. 


Part (c) above will have produced at least one more edge with 
current cost zero. Incorporate all such edges with zero cost in the 
partial graph, remove any edge which now has a non-zero cost, 
delete all labels, and return to Part A, Step 1. 
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This completes the description of the Hungarian algorithm. It is 
instructive at this point to compare this algorithm and the algorithm for 
finding a maximum matching given in the previous section with the 
maximum flow algorithm described in Networks 1. As we pointed out in 
Section 2, the problem of finding a maximum matching in a bipartite 
graph is closely related to the problem of finding a maximum flow in a 
network. The process of finding an alternating path is essentially the 
same as that of finding a flow-augmenting path in the maximum flow 
algorithm. The algorithms for these two problems find a flow or a 
matching at minimum cost, and then increase the flow or the matching in 
stages, using a procedure which ensures that the cost is a minimum at each 
stage. Such a procedure is called a least-cost flow-augmenting cycle. 


Remark Some people prefer not to draw the partial graphs when 
working through the algorithm and instead record everything on the cost 
matrices. This is very simple to do; the edges of the partial graph 
correspond to zero entries in the cost matrix, and those zero entries in the 
matrix corresponding to a matching in the partial graph are indicated by 
underlining them. The labels required during the labelling procedure can 
also be indicated on the matrix by placing them by the vertex labels. 


Use of dummy vertices 


Unlike the maximum matching algorithm and the Hungarian algorithm 
for the transportation problem, the Hungarian algorithm for the 
assignment problem as stated above works only when the two sets X and Y 
have the same number of elements. However, given a problem in which 
this condition is not satisfied, it is possible to solve the problem by 
including dummy vertices in the initial bipartite graph. The following 
worked problem illustrates the method. 


Worked problem 3.1 


A building contractor advertises five jobs — those of bricklayer, carpenter, 
decorator, electrician and plumber. There are four applicants 1, 2, 3 and 4 
— one for carpenter and decorator, one for bricklayer, carpenter and 
plumber, one for decorator, electrician and plumber, and one for carpenter 
and electrician. 


The applicants have different abilities to do the various jobs, as shown in 
the following cost matrix — the lower the cost, the more able is the 
applicant to do that particular job. 


jobs 
d 


applicants 


Use the Hungarian algorithm to assign each applicant to a job in such a 
way that the total cost is a minimum. 
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> OO 63 


This is Example 3.6 of the Introduction 
unit. 


Solution 


In order to apply the algorithm, we introduce a dummy applicant 5 with the 
same ability to do each job. Since, in the final solution, we do not wish to 
assign applicant 5 to a job, we take the corresponding cost to be larger than 
those of the genuine applicants; let us take this cost to be 20. 


jobs 
d 
1 b 
2 c 
3 d 
4 e 
oe ae ee 20 20 5 p 
cost matrix bipartite graph 
We now apply the algorithm in the usual way, ignoring the cells with no We can think of these non-entries as 
entries in the cost matrices. very large costs which do not affect the 
calculations. 


We construct the first revised cost matrix and the first partial graph. 


a of 


2 1 b 
Z 2 c 
2 3 d 
1 4 e 
ae, 0 OO O Gg. .9 5 p 


first revised cost matrix first partial graph 
By inspection, we find a matching with four edges: let us choose the edges 
id, 29, 4c, 5b. 
Part A 


We carry out the labelling procedure and do not achieve breakthrough. 


Part C 
1 b 
is) 2 e 
(x) 3 d 
4 e 
5 Pp (3) 


We revise the cost matrix as follows. 
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We find that 6 = 1. 


We obtain the following revise 


2g. 0 | 


d cost matrix and revised partial graph. 


We carry out the labelling procedure and do not achieve breakthrough. 


(d) 


(p) 


(+) 


Part C 


1 b 
2 ¢ 
3 d 
4 E 
0 1 5 p 
1 b 
2 c 
3 d (3) 
4 e 
5) -_ 


We revise the cost matrix as follows. 


We find that 6 = 1. 
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We obtain the following revised cost matrix and revised partial graph. 


1 b 
2 C 
2 d 
4 e 
5 P 


Part A 


We carry out the labelling procedure and achieve breakthrough at e. 


(dyn po 
2 c (1) 
(*) 3 d (3) 
(c) 4 = 1 
@ 4 Pp G) 


An alternating path is 3p2b5e. 


Part B 


We obtain the following revised matching. 


1 b 

y) C 

3 d 

4 € 

5 P 
Part A | 


No vertex in X can be labelled, so we STOP. 
We now ignore the dummy applicant 5. 
Thus an optimum assignment is 
1-decorator, 2—bricklayer, 3-plumber, 4—carpenter 
with total cost 
(34+0-1)+(4+0+0)+6+0-2)+(1+0+0)=10 
calculated from the final revised cost matrix, or 
2+44+3+1=10 


calculated from the original cost matrix. 


a 


3.2 Bottleneck assignment problem 


Many problems which occur in practice are variations of the basic 
assignment problem, and can be solved by modifying the Hungarian 


algorithm. We now give an example of such a problem, called the 


bottleneck assignment problem, and show how the basic algorithm can 
be adapted to solve it. 


Suppose that a product is to be made on a serial production line which 
involves a number of activities, each of which must be completed once to 
produce one item of the product. A number of people are to be assigned to 
these activities. If we know how long it takes each person to complete 
each activity, how can we assign the people to the activities so that the 
time taken to complete one item of the product is as short as possible? 


It is clear that the production line can operate only as fast as the rate of 
working of the slowest person on the line — the bottleneck. What we need 
to do, therefore, is to find an assignment of people to activities such that 
the longest time taken by any person to do his or her assigned activity is a 
minimum. In the following worked problem, we show how to solve such a 
problem by modifying the algorithm for the assignment problem so that, 
instead of finding a maximum matching at minimum total cost, we find a 
maximum matching for which the largest individual cost (in this case the 
time taken by the slowest worker) is a minimum. 


Worked problem 3.2 


Suppose that three people are to be assigned to three activities on a 
production line, and that the times (in minutes) taken for the people to do 
the various activities are given by the following cost matrix. 


activity 
WY 2 
Xy 3 2 
person % 
X3 2 2 


How can we assign-the three people to the three activities in such a way 
that the production line can operate as quickly as possible? 


Solution 


We denote the time taken to complete one item of the product by p, so that, 
in a particular assignment of people to activities, p is the longest time 
taken by any one person to do his or her assigned activity. At each stage of 
the modified algorithm for this problem, we need to keep account only of 
the parameter p. This makes the procedure considerably simpler than 
that of the algorithm for the basic assignment problem, where at each 
stage we have to keep account of the total cost, which we do by assigning 
weights to all the vertices of the graph. 


The lowest cost in the above cost matrix is 1, corresponding to worker x; 
doing activity y;. We set p initially to this value, because it cannot be 
lower, and subtract 1 from all the entries in the cost matrix. 


This gives a new cost matrix from which we can construct the 
corresponding partial graph, as we did in the assignment problem, using 
only those entries in the cost matrix which are zero. The new cost matrix 
and partial graph corresponding to p = 1 are as follows. 
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‘The time taken to complete one item 
of the product’ is ambiguous here. 
What is meant is the time interval 
between two items rolling off the 
production line. This is not the total 
sum of the times for the various 
activities, because the workers are 
working on the activities at the same 
time. 


Note that when we have a partial 
assignment, the current value of pwill 
be the longest time taken by any 
person to do the activities assigned to 
him or her so far. 


x; @ OY; 
X> @ OY? 


The partial graph has only one edge, so a maximum matching comprises 
just this edge, shown by a thick line. A complete assignment has not been 
found, so we must enlarge the partial graph. The minimum of all the non- 
zero costs is 1, so we reduce the costs of all non-zero entries in the matrix by 
1 and increase the value of p by the same amount. This gives the following 
cost matrix and partial graph, corresponding to p = 2. 


Yi Yo 3 
xy 1 0 0 xy 41 
X» 1 2 1 X72 @ Y2 
x3 0 0 0 x3 ¥3 


We can now include the edge x; y2 in the matching, so this is shown by a 
thick line. 


Again, a complete assignment has not been found, so we continue updating 
the cost matrix, finding the minimum non-zero entry (again 1) and reducing 
each non-zero entry by this amount and increasing p by the same amount. 
This gives the following cost matrix and partial graph, corresponding to 
p = 3. 


Y¥1 Yo Y3 
xy 0 0 O xy Y} 
X92 0 1 0 X2 Y2 
X3 0 0 0 x3 ¥3 


We can now include the edge x y; in the matching, so this is shown by a 
thick line. A complete assignment has now been made, and so the minimum 
time for the completion of one item is 3 minutes (the final value of p). & 


We state the algorithm formally as follows. 


Algorithm for the bottleneck assignment problem 


START witha given cost matrix. 
Set p = 0. 


STEP 1 (a) Setd=smallest non-zero entry in current cost matrix. 
(b) Reduce each non-zero entry in the current cost matrix by d. 
(c) Increase p by d. 

STEP 2 Use the matching algorithm to find a maximum matching in the 


bipartite graph. 
If a complete assignment is achieved, STOP: this is an optimum 
assignment and the value of p is the required answer. 


Otherwise, return to STEP 1. 


Note that there are other possible 
optimum matchings in the final 
partial graph corresponding to 
alternative optimum assignments. 


The simple example in the above 
worked problem is chosen to illustrate 
only the modifications to the 
assignment algorithm. Thus in Step 2 
we found a maximum matching by 
inspection each time. For a more 
complicated example, in Step 2 it is 
necessary to find an alternating path 
with respect to the existing matching 
and use it to obtain a maximum 
matching. 


of 


Problem 3.1 


Four people are to be assigned to four activities on a serial production line. 
The times (in minutes) taken for the people to do the various activities 
are as follows. 


activity 
Yo 3 
% |. 60% 4° Bee 
person *2 Pes AG tr (ge 3 
% | Bos BG ye 
Xe | Burt, ot Z 


Assign the four people to the four activities in such a way that the 
production line can operate as quickly as possible. 


3.3 Computer activities 


The computer activities for this section are described in the Computer 
Activities Booklet. 


ransportation problem 


The well-known transportation problem, sometimes called the Hitchcock 
problem after F. L. Hitchcock, who studied the problem in 1941, can be 
stated as follows. 


} 


/ 


Transportation problem 


/ 
A manufacturer has a number of factories, each of which can supply a 
particular product to a number of warehouses. Each factory has a fixed 


output, and there is a transportation cost involved in sending the 
product from each factory to each warehouse. To which warehouses 
should the products of each factory be sent so that the requirements of 
all the warehouses are satisfied at minimum total cost? 


Example 4.1 


An example of the transportation problem involving three factories and 
four warehouses is represented by the following bipartite graph. The costs 
of supplying each warehouse from each factory are given in the cost 
matrix to the right of the graph. 
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factories warehouses 


@ A 
supply (2) As | °2.: S23 --2 


Oe Mt Aa 26h Fer" 


In the bipartite graph, the factories A;, A, and A are represented by the 
black vertices, which we call supply vertices. The warehouses bj, b2, b3 
and by, are represented by the white vertices, which we call demand 
vertices. 


The circled numbers next to the supply and demand vertices on the graph, 
and also on the cost matrix, are the numbers of units of the product 
supplied by the factories or demanded at the warehouses. In this 
example, the total supply and the total demand are both 15 units. 


A flow pattern which satisfies the supply and demand constraints is 
indicated by the thick edges of the bipartite graph. The number beside 
each of these edges is the number of units of the product to be sent along 
the route represented by that edge. The other edges represent unused 
routes. We can calculate the total transportation cost for this flow pattern 
by multiplying the amount of flow along each edge by the cost associated 
with that edge. Thus the total cost is 


(4x5)+(1x«7)+(1x2)+(3x5) + (4x 4) + (2x 4) = 68. 


Although this flow pattern satisfies the flow constraints, it is not a 
minimum-cost solution. & 


The transportation problem is closely related to the assignment problem 
which we discussed in the previous section. If, in the assignment problem, 
we interpret the assignment of a worker to a job as being a ‘flow’ of one unit 
between the corresponding vertices of the bipartite graph representation, 
then we see that the assignment problem is a special case of the 
transportation problem in which the supplies and demands at the vertices 
are all 1. The transportation problem is more complicated, because we 
need to specify the amount of flow along each edge included in a solution, 
and also because a single supply vertex (factory) may send goods to more 
than one demand vertex (warehouse). 


The Hungarian algorithm for the transportation problem which we give 
in this section is based on the algorithm for the assignment problem. 


Note that for a solution to be possible, 


the total supply must be at least as 
large as the total demand. 


We ask you to find a minimum-cost 


solution by using an algorithm in 


Problem 4.1. 


4.1 Hungarian algorithm for the transportation problem 


' Before setting down the detailed procedure of the algorithm, we describe 
the method in general terms. 


The first step is to construct a partial graph, using the same procedure as 
in the algorithm for the assignment problem. We then assign the 
maximum possible flow to the edges of the partial graph, subject to the 
supply and demand constraints at each vertex. The way in which the 
partial graph is constructed ensures that this is a minimum-cost flow — in 
other words, we cannot obtain the same value of the total flow at a lower 
cost by using edges not in the partial graph. If the flow thus obtained does 
not satisfy the demands at all the demand vertices, then we enlarge the 
partial graph, using the same procedure as in the algorithm for the 
assignment problem. We can then increase the value of the total flow by 
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using the new edges. We continue in this way until we eventually obtain a 
minimum cost solution in which all the demands are satisfied. It is 
important to note that, as with the assignment problem, there may be 
other minimum-cost solutions. 


The Hungarian algorithm thus builds up the flow in stages, ensuring that 
at each stage the total flow is at a minimum cost — thus it uses a least- 
cost flow-augmenting cycle. 


In the Hungarian algorithm for the transportation problem, we are not 
seeking a matching, as we were in the algorithm for the assignment 
problem. (This is because there may be flows from one supply vertex to a 
number of demand vertices.) So, instead of looking for alternating paths, 
we look for paths in which there is a flow on even edges — that is, on the 
2nd, 4th, 6th, ... edges. There may or may not be a flow on odd edges of the 
path. 


In the following audio-tape session, we illustrate the use of this 
algorithm by working through an example. Before listening to the tape, 
read through the description given below. Have this description and the 
Audio-tape Notes to hand when you listen to the tape. 


Hungarian algorithm for the transportation problem 


START with no flow in a bipartite graph whose vertices are divided 
into two sets: the supply vertices A,,..., A, and the demand 
vertices D1, ..., Dy. 


Construction of initial partial graph 


STEP 0a Assign weights to the supply and demand vertices as follows. 


To each supply vertex, assign a weight equal to the lowest cost on 
any edge incident with that vertex. 


Decrease the cost on each edge by the weight on the supply 
vertex with which it is incident. 


Repeat this for each demand vertex. 


STEP Ob Construct a partial graph consisting of the vertices of the 
original bipartite graph, together with only those edges whose 
current cost is Zero. 


Part A: labelling procedure 
STEP 1 Label with (*) each supply vertex whose supply has not all been 
allocated to edges. 


If no such vertex exists, STOP: the current solution is a minimum- 
cost solution. 


Otherwise, go to Step 2. y 


STEP 2 Choose a newly labelled supply vertex Aj, and label with (Aj) 
all unlabelled demand vertices which are joined by an edge of 
the partial graph. 


Repeat this for all newly labelled supply vertices, and then go 
to Step 3. 


STEP 3 Choose a newly labelled demand vertex b;, and label with (b/) 
all unlabelled supply vertices which are joined to b; by an edge 
which has been allocated a flow. 


Repeat this for all newly labelled demand vertices. 
Repeat Steps 2 and 3 until 


EITHER a demand vertex whose demand is not satisfied is labelled (this 
is called breakthrough), in which case go to Part B, Step 4, 


OR no more labelling is possible, in which case go to Part C, Step 6. 
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Such a path is analogous to a flow- 
augmenting path used in the 
maximum flow algorithm in 
Networks 1. 


Part B: flow-augmenting procedure 


STEP 4 Start with the breakthrough vertex, and trace back through the 

labels until a supply vertex which is labelled (*) is reached. A 

flow-augmenting path has been found. The flow-augmenting path plays a 
similar role in this algorithm to that of 


Calculate the maximum flow F which can be sent along this path 4. sfetnduine oath i the alvorittim for 


as follows. the assignment problem. 

First, find the minimum of all flows assigned to EVEN edges of the 

path (that is, the 2nd, 4th, ...edges). This is the maximum 

possible backward flow. 

Next, find the minimum value of: Note that for a flow-augmenting path 


with just one edge, there is no 
backward flow, so we just take the 
e the required demand at the end vertex; minimum of the supply and the demand. 


e the available supply at the start vertex of the path; 


e the maximum backward flow. 
This is the required value of F. 
STEP 5 Forma new flow pattern as follows. 


(i) Increase the flows on the ODD edges of the path by 
this value of flow, F. 


Decrease the flows on the EVEN edges by F. 

(ii) Decrease the available supply at the START vertex by F. 
Decrease the demand at the END vertex by F. 

Remove all the labels and return to Part A, Step 1. 


Part C: modification of partial graph procedure 
STEP 6 Find the lowest current cost on any edge of the original bipartite 
graph which 
e starts at a LABELLED supply vertex; 
e ends at an UNLABELLED demand vertex. 
Call this lowest cost 6. 
STEP 7 (a) Increase the weight on each labelled supply vertex by 6. 
(b) Decrease the weight on each labelled demand vertex by 6. 


(c) For each edge which joins a labelled supply vertex to an 
unlabelled demand vertex, decrease the current cost by 6. 


(d) For each edge which joins an unlabelled supply vertex to 
a labelled demand vertex, increase the current cost by 6. 


Incorporate all edges which now have a current cost of zero in 
the partial graph, remove any edge which now has a non-zero 
cost, delete all labels, and return to Part A, Step 1. 


€> Go €3 
ee 
Problem 4.1 
Use the Hungarian algorithm to find a minimum-cost solution to the 
transportation problem described in the text at the beginning of this ®© A, 
section, and calculate the total transportation cost. The cost matrix is 
repeated in the margin. depp EAP URN SG SS Fey 


OF a ae Pe oe 
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Use of a dummy vertex 


In each of the examples of the transportation problems which we have 
considered, the total supply is equal to the total demand. In other cases, 
the total supply may be greater than the total demand — for example, in 
the previous problem it may be a simple matter to increase the supply of 
each factory by 2 units. To enable us to use the Hungarian algorithm to 
solve such problems, we introduce a dummy vertex. This represents an 
artificial warehouse (or other demand vertex) whose demand is chosen to 
be equal to the excess supply. Thus, for the example of Problem 4.1, if each 
factory can supply 2 additional units and the warehouse demands are the 
same, we represent the problem by the following bipartite graph and cost 
matrix. 


supply demand 


@ Ai b, 6) ® A 


(6) Ag bo(3) sapply @a,/ 2 °5 3 #2 0 


(8) As b3(5) it G&G. 3 4 ae 


bs(6) 


Vertex bs is the dummy vertex with an artificial demand of 6, equal to the 
excess supply. In the cost matrix, we must set the transportation costs from 
each supply vertex to this dummy vertex to be equal — we have chosen 
here to set them all equal to 0. We can then apply the algorithm in the 
usual way. The flow pattern produced will send 6 units to the dummy 
vertex. To obtain a flow pattern which is a solution to the original 
problem, we simply remove all flows to the dummy vertex. 


4.2 Transhipment problem 


In the transportation problem discussed above, we assumed that the 
products from each factory are sent directly to one or more warehouses. 
However, in a practical situation, it may be cheaper to send the products 
from a factory first to an intermediate point (one of the other factories or 
warehouses) and then on to the destination warehouse. 


The sending of goods via intermediate points is called transhipment, and 
a problem in which transhipment is allowed is called a transportation 
problem with transhipment, or, more simply, a transhipment problem. The 
graphical representation of a transhipment problem is not a bipartite 
graph, since we can send goods from factory to factory, or from warehouse 
to warehouse. However, we can easily transform the graphical 
representation into a bipartite graph — to do this, we represent each 
factory or warehouse by two vertices, one representing the factory or 
warehouse as a supply point, and one representing it as a demand point. To 
illustrate how this works, we consider a specific example. 


Example 4.2 


Suppose we have a transportation problem for which the transportation 
costs are as shown in the following cost matrix, in which the circled 
numbers are the amounts of supply or demand at the vertices. 
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@) A, 
supply ©) A,| 17 14 11 


@) Aa 18 18" 15 


If transhipment is allowed, then there are costs associated with the 
various transhipment routes. We suppose that these are as given in the 
following cost matrix. 

demand 


® @® ® ®© © 


Aa, &y b, 


@ A, 
© A, 
@ As 
@ B, 
@® By 
© Bs 


In this matrix, A; represents the first factory considered as a supply 
vertex, and a, represents the same factory considered as a demand vertex. 
Similarly, B, represents the first warehouse considered as a supply 
vertex, and b, represents that warehouse considered as a demand vertex. 
Since we must consider each factory and warehouse as both demand and 
supply vertices, we have a total of six supply vertices and six demand 
vertices. Note that the part of the cost matrix for the original vertices 
Aj, Az, A3 and b,, bz, bz is the same as the original cost matrix. 


supply 


12. 17 0 2 


i -4a 3S 4 2 0 


As before, the circled numbers indicate the amounts demanded or 
available at the vertices considered as demand or supply vertices. 
Factories considered as demand vertices have zero demand, and 
warehouses considered as supply vertices have zero supply. The fact that 
the cost matrix is not symmetric means that outgoing shipments and 
returning shipments do not cost the same. The entries on the diagonal of 
the matrix are all zero, since they are the costs of transporting goods 
between supply and demand vertices which represent the same factory or 
warehouse. £ 


By splitting the vertex representing each factory or warehouse into two 
vertices — a supply vertex and a demand vertex — we have cast the 
transportation problem with transhipment into the same form as a 
transportation problem without transhipment. We do this so that we can 
apply the Hungarian algorithm for the transportation problem to the 
problem with transhipment. However, before we can apply the 
algorithm, a further modification is required. The reason for this is the 
following. When we apply the Hungarian algorithm, we can ship goods 
from a demand vertex to a supply vertex only if there is already a flow 
from that supply vertex to that demand vertex. This means that, for the 
transportation problem with transhipment, we cannot ship goods from a 
particular factory considered as a demand vertex to the same factory 
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considered as a supply vertex. To overcome this problem, we add to all the 
supplies and demands an amount of demand or supply at least as great as 
the total demand or supply. This ensures that any vertex considered as a 
demand vertex can accommodate any transhipment which could occur in 
the solution to the problem. When the algorithm has been completed, we 
remove these artificial supplies and demands (and the corresponding 
artificial flows) to obtain a solution to the original problem. 


Worked problem 


Use the Hungarian algorithm for the transportation problem to obtain a 
minimum-cost solution to the above transhipment problem. 


Solution 


The total supply and the total demand are both equal to 15. We add an 
amount greater than this to all the supplies and demands — we choose 20 
as the amount to be added. 


When we apply the algorithm, we form a partial graph, using the cost 
matrix given in the above text. Since all the entries on the main diagonal 
are zero, and all other entries are non-zero, this partial graph consists 
only of the edges joining A; to a,, Az to a2, and so on, as follows. 


Q2) Ay a1 Q0) 
26) Az a2 Q0) 
7) As a3 Q0) 
Q0) By b1 Q5) 
20) Be bo 24) 
20) Bs bs Q6) 


The labelling procedure and the procedure for modifying the flow pattern 
result in a flow of 20 being sent along each edge — these are the artificial 
flows necessary to ensure that any transhipment can be accommodated. 
The resulting flow pattern is shown below. The circled numbers next to the 
vertices represent unallocated supplies or unsatisfied demands — these 
are just the original supplies and demands. 


QA a4) 

6) Ay 20 a>(0) 

DAs $§oa@ 

Oe anne 

()B, 20 b,(4) 

@B,-# bs (6) 
We continue to apply the algorithm in the same way as for a 
transportation problem without transhipment: the procedure is exactly 


the same as that used in the first worked problem. We thus obtain the 
following flow pattern. 
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©6666 
cea 6 


5 b3(0) 


The flow path A3a,A,b, involves transhipment via the factory A,. To 


accommodate the flow of 2 along this path, the artificial flow along the 
edge Aj,a, has been reduced from 20 to 18. No other flow path involves 


transhipment. 


To obtain the flow pattern which is a solution to the original problem, we 
must subtract all the artificial flows of 20 from the flow pattern given 
above. This gives the following solution, in which we have indicated the 
direction of flow along each edge by an arrow. The original supplies and 
demands are shown next to the vertices. 


(2) Ay . 


4.3 Computer activities 


The computer activities for this section are described in the Computer 
Activities Booklet. 


Further reading 


Discussion of the material in Sections 1 and 2 is included in many books on 
combinatorics — for example: 


J. O. Clark and D. A. Holton, A First Look at Graph Theory, World 
Scientific Publishing, 1991; 


G. Chartrand and O. R. Oellerman, Applied and Algorithmic Graph 
Theory, McGraw-Hill, 1993; 


R. J. Wilson, Introduction to Graph Theory, 4th edition, Longman Group, 
1996; 


C. L. Liu, Introduction to Combinatorial Mathematics, McGraw-Hill, 1968. 


Most books on network theory include a section dealing with aspects of the 
topics of Sections 3 and 4 — for example: 


L. R. Ford, Jr. and D. R. Fulkerson, Flows in Networks, Princeton University 
Press, 1962; 


E. Lawler, Combinatorial Optimization: Networks and Matroids, Holt, 
Rinehart and Winston, 1976. 
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Exercises 


Section 1 


Matching problem 


1.1. Five people have each won a holiday to one of five resorts — 
Blackpool, Cannes, Edinburgh, Paris and San Francisco. Of the winners, 
one will go only to Britain, one will not go to France or Scotland, two will 
not go to Britain, and the fifth will go only to San Francisco. 


(a) Draw the bipartite graph representing this situation. 

(b) By finding a maximum matching in this graph (by inspection), 
determine whether it is possible for each winner to go to a different 
resort. 


1.2 A computer agency advertises for a systems analyst (S), a 
programmer (P) and a computer manager (M). There are four applicants — 
one is a qualified systems analyst (A,), two are both qualified systems 
analysts and programmers (A, and A;3), and the other is a qualified 
programmer and computer manager (A,). 


(a) Draw up a table showing, for each possible subset of jobs, the number 
of jobs in that subset and the number of suitable applicants for that 
subset of jobs. 


(b) For each subset, indicate whether the marriage condition is 
satisfied, and decide whether it is possible for all three jobs to be 
filled. 


(c) Draw the bipartite graph representing this situation and find, by 
inspection, a maximum matching in this graph corresponding to the 
maximum number of jobs that can be filled. 


1.3 Does there exist a maximum matching containing 4 edges in the 
following bipartite graph? If so, find such a matching; if not, explain why 
such a matching does not exist. 


*y Yi 
x2 ¥2 
“5 ¥3 
X4 Y4 

Ys 


1.4 Five tutors T,,T2, 73,14, T5; apply to teach five courses; just one tutor 
is required for each course. To find out how many of the tutors can be 
appointed to courses which they are competent to teach, the marriage 
theorem is applied to each subset of tutors. It is found that the marriage 
condition is satisfied for all except the following subsets: 

tutors T,, T, and T; are between them competent to teach only 

2 courses; 
tutors T; and T; are between them competent to teach only 1 course. 


(a) How many of the five tutors can simultaneously be appointed to 
courses which they are competent to teach? Explain your answer. 


(b) For which of the following subsets can all the tutors in the subset be 
appointed to courses which they are competent to teach? Explain 
your answer. 


wv 


(1) {T,,T2,T3, T;} i tT} 
(3) {T;, Tp, T4, Ts} / (4) i; T3, Ts} * 
(5) (T1, To, T3, Ts} és (6) te T3, ta; Ts} ¥ 


ay 


1.5 A university has 5 vacancies V,,V 2,V3,V4,V5, and there are 
5 applicants A,, Az, A3, Ag, As, whose suitability for the different jobs is 
shown in the following table. 


job -—_ suitable applicants — 
ae oe 
V2 A3, Ag 

V3 A, 

Vi — 

V5 Ay, Az, As 


(a) Draw a bipartite graph to represent this situation. 


(b) List all the subsets of vacancies for which the marriage condition is 
not satisfied. 


(c) For each of these subsets, find the value of the expression in the 
modified marriage theorem, and hence find the largest number of 
vacancies which can be filled. 


(d) Find, by inspection, a matching in the bipartite graph with this 
number of edges. 


Section 2 


Maximum matching algorithm 


2.1 Find, by inspection, an alternating path in the following bipartite 
graph, and hence obtain a new matching with one more edge. 


x4 Y 
X> Y2 
X3 ¥3 
X4 Y4 


2.2 Carry out the labelling procedure for the graph in Exercise 2.1. 
2.3 Consider the following labelled bipartite graph. 


a | Y1 
(y4) Xp Yo (x4) 
(y3) X3 ¥3 (Xp) 
(Ys) x4 Y4 (x4) 
(+) x5 @—— Ys (Xs) 


Starting from the breakthrough vertex, trace back through the labels to 
find an alternating path, and hence obtain a new matching with one more 
edge. 


2.4 Use the maximum matching algorithm to find a maximum matching 
in the following bipartite graph, starting from the given initial 
matching. 


x4 Y1 
x2 Y2 
X3 ¥3 
X4 Y4 
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Section 3 


Hungarian algorithm for the assignment problem 


3.1 Construct the first revised cost matrix from the original cost matrix 
given below, and draw the corresponding partial graph. 


3.2 For the following partial graph and current cost matrix, find the 
lowest non-zero current cost 6 and construct the revised cost matrix. 


i oe ee f= 2 


ao a: a a | x4 ¥1 
o_o a a Xp Yo (x4) 
4%; 3 Q_2 4 (Y2) x3 Y3 

—2—%4-+_5——_)—_ 4 () 4 Y4 


3.3 Four candidates apply for four jobs. The selection committee rates 
the suitability of each candidate for each iob and expresses this as a cost, 
where the lower the cost, the more suitable is that candidate for that job. 
Allocate the candidates to the jobs in such a way that the total cost is as 
small as possible, using the Hungarian algorithm. The costs are as 
follows. 


jobs 
ln ye ye Ms 
% | Sto Io _ 84 
y-| Ts 4&2 2 
applicants 
x3 8 8 5 
X4 8 8 5 


Bottleneck assignment problem 


3.4 Four people are to be assigned to four activities on a serial production 
line. The times (in minutes) taken for the people to do the various jobs are 
as follows. 


activities 
YM Wo Ys Yo 
xX 1 3 4 8 
X | 4 4 => 2 
people : 7 3 — 
X4 6 4 be 3 


Use an appropriate algorithm to assign the four people to the four 
activities in such a way that the production line can operate as quickly as 
possible. 


ao 


Section 4 


Hungarian algorithm for the transportation problem 


4.1 Given the following partial graph of a transportation problem, 
carry out the labelling procedure. 


2pb1(0) 


4.2 Three factories, A;, Az and A3, can supply 10, 15 and 20 units of 
product, respectively. The demands for these products at four warehouses, 
b,, b>, bz and by, are 5, 12, 13 and 15 units, respectively. The costs of 
transporting one unit of product from each factory to each warehouse are 
given in the following cost matrix. 


Find a minimum-cost solution to this transportation problem, using the 
Hungarian algorithm. 


4.3 Describe briefly how you would modify your solution to Exercise 4.2 
if the supply at factory A, were 18, rather than 15. 


Transhipment problem 
4.4 A transhipment problem has the following cost matrix. 


demand 


® © ® © 


Qa, a, dy bo 


@) A, 

(3) Ay 
supply 

@) By 


©) By 


Use the Hungarian algorithm to find a minimum-cost solution, and 
determine the value of this minimum cost. 


40 


Solutions to the exercises 


1.1 
(a) The bipartite graph representing this situation is: 


P, Blackpool 
P» Cannes 

P3 Edinburgh 
Fy Paris 

Ps San Francisco 


(b) There are only two maximum matchings: 


Py B P, B 
P, C P, C 
P; E P, E 
Cpe P, P 
P.o———-__——0 S P,o————05 


These matchings show two ways in which all five winners can go to 
different resorts. 


+2 


(a) 
number number of 


subset of jobs applicants of jobs applicants 
{S, P, M} A;, A, A3, Ag 3 + 
iS, Fi A;, Az, A3, Ag z 4 
{S, M} Ay, Az, A3, Ag 2 “ 
{P, M} Ao, A3, Ag 2 3 
{S} A;, Ao, Ag 1 3 
{P} Az, Ag, Ag 1 ’ 
{M} Ag 1 1 


(b) The marriage condition is satisfied for every subset, so, by the 
marriage theorem, all three jobs can be filled by suitable applicants. 


(c) The bipartite graph representing this situation is shown below. A 
maximum matching of three edges is indicated by three thick lines. 


Ay 

s 
Ay 

P 
Ag 

M 
Ag 


1.3 There is no maximum matching containing four edges in this bipartite 
graph because each of the vertices x,, x3 and x4 is joined only to the vertices 


y3 and ys. 


1.4 


(a) Only four tutors can simultaneously be appointed to courses which 
they are competent to teach, because only T, Ty, T; and one of T; and 
T3 can be appointed. 


(b) The sets containing only one of T, and T; — namely, sets (2) and (3). 


Other matchings with three edges are 


possible. 
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hie 
(a) 


(b) 


(c) 


(d) 


The bipartite graph representing this situation is: 


V; Ay 
V> a 
V3 A3 
V4 A, 
: Ae 


The subsets of vacancies for which the marriage condition is not 
satisfied are listed below: 


number of number of 
vacancies applicants vacancies m applicants p 
{V2, V3, Va} As, Ag 3 2 
{V3, V4} Ag Z 1 


For each of these subsets, the expression p + (n — m) in the modified 
marriage theorem has value 4, and so the largest number of 
vacancies that can be filled by suitable applicants is 4. 


Both subsets of vacancies for which the marriage condition is not 
satisfied contain the vacancies V3 and V,. Therefore, if we remove 


one of the vacancies V3 and V4, the remaining vacancies can be filled. 


A solution is shown in the following diagram. 


yf =O A, 
V2 OAs 
ves Ay 
y ———oO A, 


V;o———_—0 A, 


2.1 The only alternating path is x3 y; x1 Y3 X4 Y4. 


To obtain the new matching, we take the edges of the alternating path 
NOTIN M — that is, x3 y;, x; y3 and x4 yg — together with the edges of M 
NOT IN the alternating path — that is, x, y,. This gives the matching 


shown in the margin. 


2.2 Wecarry out the following steps. 


STEP 1 


Siar 2 


STEP 3 


Sine 2 


SrEP 3 


STEP 2 


Label with (*) each vertex in X which is not incident with any 
edge in the current matching M — that is, x3. 


Label with (x3) all unlabelled vertices in Y joined to x3 by an edge 
NOT IN M — that is, y4. 

Label with (y;) all unlabelled vertices in X joined to y, by an edge 
IN M — that is, x. 

Label with (x;) all unlabelled vertices in Y joined to x; by an edge 
NOT IN M — that is, ys. 

Label with (y3) all unlabelled vertices in X joined to y3 by an edge 
IN M — that is, %. 


Label with (x4) all unlabelled vertices in Y joined to x, by an edge 
NOT IN M — that is, y. and y4. | 


Breakthrough occurs at y4, so the labelling procedure terminates. 
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Other solutions are possible. 


| Y} 

Xo Y2 

- Y3 

X4 Ya 
(yy) xy Y, (x3) 

x2 Yo (X4) 
(+) x3 43 (x) 
(ys) Xq Yq (X4) 


Note that if we trace back through the 
labels from the breakthrough vertex 
y4 we obtain the alternating path in 


Solution 2.1. 


2.3 The breakthrough vertex is y2. Tracing back through the labels from 
Y2, we obtain the alternating path x5 ys X4 Yo. 


The new matching is obtained from the existing one by replacing the edge 


X4 Ys by the edges x4 y. and xs ys, as follows. 


2.4 Applying the labelling procedure, we obtain one of the following 


labelled graphs. 


(yp) % y 

x2 Y2 
(+) X3 ¥3 
(+) %4 Y4 


(a) 


In each case, breakthrough occurs at y3. 


x4 
x2 
X3 
X4 


x5 


(x3) 
(x,) 


(x1) 


yy 


Y2 


¥3 


Y4 


Ys 


(y)) 


(+) 


(+) 


In case (a), an alternating path is x3 y; x1 y3. 


In case (b), an alternating path is x4 y; x; y3. 


xy yy, (4) 
x2 Y2 (x4) 
x3 a WYg. (xy) 
X4 ¥4 

(b) 


We thus obtain one of the following new matchings. 


x4 Y1 
x2 Y¥2 
X3 Y¥3 
X4 Y4 


(a) 


Applying the labelling procedure again, we obtain one of the following 


labelled graphs. 


x} Yj 

x2 Y2 

(yy) x3 ¥3 
(+) x4 Y4 


(a) 


(x4) 


(+) 


(y,) 


xj Y 

x2 Y2 

x3 ¥3 

X4 Ya 
(b) 


xy y, (%3) 
x2 Y2 
X3 Y¥3 
X4 Ya 


(b) 


Breakthrough is not achieved in either case, so we STOP. 


Thus a maximum matching has three edges, as shown above. 


The point of the exercise is to 
demonstrate how the algorithm 
terminates. However, in this small 
example it is clear from the graph 
that a maximum matching has only 
three edges, since the vertices A3 and 
Ay, are joined only to by. 
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3, SO We give xX, weight 
3 and reduce the cost of all the entries in row x, by 3. Similarly, we reduce 
LOWS Xz, x3 and x, by their lowest costs 1, 2 and 1, respectively. The matrix 


becomes 


Next we consider the columns. Column y, has lowest cost 1, column y, has 
lowest cost 3, and columns y3 and y, each have lowest cost 0. For each 
column, we subtract the lowest cost from each entry in that column, and the 
revised cost matrix becomes 


2 xy Y1 
1 x2 Y2 
2 x3 ¥3 
1 % 1 o* 3 0 X4 Ya 
first revised cost matrix first partial graph 


The corresponding partial graph, consisting of all the vertices together 
with the edges with zero current cost, is shown above on the right. 


3.2 We start with the following situation, and apply Part C of the 
Hungarian algorithm for the assignment problem. 


a oo: 0 


5 x4 Y1 
3 x2 Yo (4) 
4 (yo) X3 Y3 
2 X4 Bit tc Orit 5 (*) x4 Y4 


STEP 6 We draw a horizontal line through each labelled vertex in X. 
We draw a vertical line through each labelled vertex in Y. 


The smallest entry with only a horizontal line through it is 
min G, 2,4, 5),so6=2. 

STEP 7 We decrease all entries with only a horizontal line through 
them by 2, and increase the weights on the corresponding 
vertices in X by 2. 


We increase all entries with only a vertical line through them 
by 2, and decrease the weights on the corresponding vertices in Y 
by 2. 
We thus obtain the revised cost matrix. 

0 -1 2 90 


3.3. We apply the algorithm as follows. 
We construct the first revised cost matrix and the first partial graph. 
- 8 - 8 9a 


Y¥3 
1%) @ 6 4 9 x v1 
oe ee % 7 
2 ae ee x, a = 
2%/14 6 3 O x4 Re 


Part A 


By inspection, we find a matching with three edges: let us choose x, y1, 
Xo Y3 and x3 y4. We label the first partial graph as follows. 


| eee Y} 

x2 Y2 

(ys) Xs ie Y3 
(%) x4 fies. Y4 (X4) 


Breakthrough is not achieved, so we go to Part C. 


Part C 
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We obtain the following revised cost matrix and revised pa 


2 i 
Ya 
1 x, 0 GO #4 -T2 


2 X2 a 4 


5 X3 i-* 3. ee 
5X4 i ee ee 
Part A 
We label the revised partial graph as follows. 
(Ys) x2 Y2 
(y4) X3 Y3 (x4) 
(+) X4 Y4 (X4) 


Breakthrough is not achieved, so we go to Part C. 
Part C 


2) Y2 
x3 ¥3 


X4 Y4 


We find that 6 = min (3, 2, 1) = 1. 


We obtain the following revised cost matrix and revised partial graph. 


Y3 Ya 


6 X4 GC 2 ee 
Part A 
We label the revised partial graph as follows. 
(y) x, Y, (x4) 
(y3) x2 Yo (x3) 
(Ys) x3 y3 (%4) 
(+) x4 Y4 (X4) 
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xy yi 
x2 Y¥2 
ss Y¥3 
X4 ¥4 


We have breakthrough at vertex yp. 


Tracing back through the labels from y2, we obtain the alternating path 
X4 V1 %4 Y2. 


We form a new matching as follows. 


1 Yi 
a Y2 
- Y¥3 
X4 Y4 
A complete assignment has been found at a total cost of Other matchings with this lowest 


value for the total cost are possible. 
14+242482 13, 


obtained from the original cost matrix, or 
(1+0+0)+(G+0-1)+(+0-4)+(6+0+ 2) =13, 


obtained from the final revised cost matrix. 


3.4 The lowest cost in the cost matrix is 1, so we set p = 1 and subtract 1 
from all the entries in the cost matrix. This gives the following cost 
matrix and partial graph, corresponding to p = 1. 


Xy 0 Z 3 7 X O——$——_— 1 


xX 3 2 2 2 X7@ OY? 


x3 6 3 2 0 X3 0 ¥3 

1 s- 3 2 = solange 
The partial graph has only two edges, joining different pairs of vertices, 
SO a maximum matching comprises these two edges, shown by thick lines. 


A complete assignment has not been found, so we repeat the process. The 
lowest non-zero cost is 2, so we subtract 2 from all the non-zero entries in 
the cost matrix, and increase p by 2. This gives the following cost matrix 


“» 


and partial graph, corresponding to p = 3 


X4 3 1 2 0 x4 Y4 


An alternating path is x4 y4%3y3, which leads to the matching 


X1Y1, X3Y3, X4Ya. 


A complete assignment has not been found, so we repeat the process. The 
lowest non-zero cost is 1, so we subtract 1 from all the non-zero entries in 
the cost matrix, and increase p by 1. This gives the following cost matrix 
and partial graph, corresponding to p = 4. 
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xy Yi 


XQ Y2 
x3 ¥3 
X4 2 0 1 0 x4 4 


An alternating path is x2 y; x; y2, which leads to the complete assignment 


XyYo, XQYy, X33, %XqYg- 


The minimum time for completion of an item is 4 minutes (the final value 
of p). 


4.1 We start with the following situation. 


2Pb1@) 


We carry out the following steps. 
STEP 1 The supplies at the supply vertices A, and A3 have not all been 
allocated, so we label A; and A3 with (+). 


STEP 2. The newly labelled supply vertex A, is joined to b, and the newly 
labelled supply vertex A3 is joined to b3, so we label by with (A;) 
and b3; with (A3). 


STEP 3. The newly labelled demand vertex b, is joined to the unlabelled 
supply vertex A, by an edge which has been allocated a flow, so 
we label A, with (b2). 


All supply vertices are now labelled. 


STEP 2. The newly labelled supply vertex Az is joined to the unlabelled 
demand vertices b,, b,and b;. We label b,; with (Az) and by, 
with (A>). The demand at b, is not satisfied, so breakthrough 
occurs at b, and the labelling procedure terminates. 


bi(0) (Ag) 
(*) Q)A bo(0) (Ay) 
(b2) (0) A2 b3(0) (As) 
(*) ()A; ba(S) (Ad) 

bs(3) 
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4.2 We start with the following situation. 


demand 
mf me nd - supply demand 
1 2 3 4 1 
Ais w- 2 30 Ay 
by (12) 
supply 4) A,|} 5 15 20 10 (15) Aa 
b3 (13) 
Q@) A;! 80 6 O 2 20) As 


bg (15) 


We construct the first revised cost matrix and first partial graph. 


Go 66h. 8 2 


biG) 


-_ (0) A1 
by (12) 
5 A> (5) Ao 
b3 (13) 
0 A3 20) A3 
by (15) 
Parts A and B 
Following the procedure of Steps 1-5, we find a maximum flow in the 
above partial graph. We thus obtain the following flow pattern. Other maximum flows are possible. 


The edges carrying a flow are shown as thick lines. We have indicated 
the amounts of flow, and have subtracted these amounts from the supplies 
and demands at the appropriate vertices. 


b(0) 
(0) Ay 10 


5 by (2) 
QQ) Ar 


Part A 
We label A, with (*) and b, with (A>). 
The partial graph with these labels is shown below. 


bi(0) (Ag) 
bp(2) 
b3(Q) 
ba(8) 


No more labelling is possible, and breakthrough has not been achieved, so 
we go to Part C. 
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Part C 
We find that 6 = min (2, 15, 3) = 2. 
We thus obtain the following revised cost matrix and partial graph. 


—2 8 QO 2 


b; 
(0) A; a0 
5 b (2) 
0) Ap 
13 bs@ 
(0) 4s 
yf 
b,(8) 
Part A 
We label the partial graph as follows. 
b,(0) (Ay) 
(0) A,el0 
: by(2) (Ag) 
(+) (10) A2 
13 BO) 
@)As 
z 
b,(8) 
Breakthrough is achieved at bz, s0 we go to Part B. 


Part B 
A flow-augmenting path is A,b,. There is no flow along this edge. 


The available supply at the start vertex A, is 10 and the required demand 
at the end vertex b, is 2, so we send a flow of 2 along A2bp. 


This results in the following flow pattern. 


b,(0) 
(0) A; a0 
5 b»(0) 
@OAKZ 
13 bs) 
@As 
7 
b,(8) 
Part A 
We label the partial graph as follows. 
bi(0) (Ag) 
(bo) (0) Aeal0 
5 b»(0) (Ap) 
(+) (8) Ap 5) 
= b3(0) (Ay) 
(b3) (0) A3 
e b4(8) (A3) 
Breakthrough is achieved at b4, so we go to Part B. 


Part B 


Tracing back through the labels from b4,, we find the flow-augmenting 
path A, b, A,b; A; dq. 
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The even edges are A, b, and A3b3 with flows of 10 and 13 respectively, 
and min (10, 13) = 10. Thus the maximum possible backward flow is 10. 


The available supply at Az is 8. 
The required demand at by, is 8. 
The maximum flow F is therefore equal to min (10, 8, 8) = 8. 


Sending a flow of 8 along the flow-augmenting path A,b, A,b3A3b4,, we 
obtain the following flow pattern. 


b, (5) 


4.3 We introduce a dummy vertex b; with demand equal to 3, the excess 
supply. The vertex bs; is joined, in the bipartite graph, to the supply 
vertices A,, Az, and A3, and the corresponding entries in the cost matrix 
are set equal to each other. We then apply the algorithm in the usual 
way, and obtain the solution of the original problem by simply ignoring 
all flows to the dummy vertex. 


4.4 We start with the following situation. 


® © ® ®@ 
@) A, 


omar oe  £°°6 


supply 
(0) B, 6 4 0 1 


on? 5. 2-4. 0 


The total supply and the total demand are both equal to 5. We add 5 to all 
the supplies and demands. 


When we apply the algorithm, we form a partial graph, using the cost 
matrix given above. Since all the entries on the main diagonal are zero, 
and all other entries are non-zero, this partial graph consists only of the 
edges joining A, to a,, Az to a, and so on, as follows. (We have omitted the 
circled numbers representing supplies and demands from the cost matrices.) 


e-- 2 -2 8 


a8) 
a(S) 
8) 
bo 7) 


Dar 
@)A2 
©Bi 
©)B. 


51 


Parts A and B 


This results in a flow of 5 being sent along each edge — these are the 
artificial flows necessary to ensure that any transhipment can be 
accommodated. The resulting flow pattern is shown below. The circled 
numbers next to the vertices represent unallocated supplies or unsatisfied 


demands — these are just the original supplies and demands. 


Part A 


on 
OA 
Os, 
Os, 


We label A; and A), with (+). 


We label a, with (A,) and a, with (A>). 


No more labelling is possible, and breakthrough has not been achieved, so 


we go to Part C. 
Part C 


3 


5 


2 


5 


a) 
a(0) 
b@) 
bo(2) 


We find that 6 = min (3, 4, 6) =3. 


We thus obtain the following revised cost matrix and partial graph. 


Part A 


—3 


We label A, and A, with (*). 
We label a; and b, with (A;). 


(2) Ay 2 a,(0) 


Breakthrough is achieved at b;,s0 we go to Part B. 


Part B 


We send a flow of 2 along the flow-augmenting path A; by. 


This results in the following flow pattern. 


ays 


Part A 
We label A> with (*). 


We label a, with (A>). 
No more labelling is possible, and breakthrough has not been achieved, so 


we go to Part C. 


Part C 


We find that 6 = min (1, 3) = 1. 
We thus obtain the following revised cost matrix and partial graph. 


a 


oe 
4 Agel 0 
ee 
0 Bee 


Part A 
We label A, with (*). 


We label Ay, a2 and b; with (A>) . 
Breakthrough is achieved at b,, so we go to Part B. 


Part B 


We send a flow of 1 along the flow-augmenting path A2),. 


This results in the following flow pattern. 


Part A (0) A; 5 a,(0) 
We label A> with (+). 


We landie ae (A). @ Ar 2042) 


We label A, with (a,) and B, with (b,). 


No more labelling is possible, and breakthrough has not been achieved, so 
we go to Part C. 


(0) Bje—= b,@) 
@B,-—~ b (2) 


Part C 


0 By 
We find that 6 = min (1, 2) = 1. 


We thus obtain the following revised cost matrix and partial graph. 


4 +5 -l 0 


we a et ee OA 1) 
5 A| 0 Gea (2), 202, (0) 
1 
i B, | oe See © Bie ©) 
0 B,.) & 2 ew (0) B,e— bo (2) 
Part A 


We label A> with (*). 

We label a,, a, and b, with (A>). 

We label A; with (a,) and B, with (b;). 
We label by with (A). 


Breakthrough is achieved at b2,s0 we go to Part B. 


Part B 


Tracing back through the labels from bj, we find a flow-augmenting There are other flow-augmenting 
path Ay a, A, by. paths. 


The even edge A, a, has a flow of 5. Thus the maximum possible backward 
flow is 5. 


The available supply at Az is 2. 

The required demand at b, is 2. 

The maximum flow F is therefore min (5, 2, 2) = 2. 

We send a flow of 2 along each of the odd edges Aja, and A;bp. 
We decrease the flow on the even edge Aa; from 5 to 3. 


We thus obtain the following flow pattern. 
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The flow path A24a,A,b, involves transhipment via the factory A,. To 


accommodate the flow of 2 along this path, the artificial flow along the 
edge A,a, has been reduced from 5 to 3. No other flow path involves 


transhipment. 


To obtain the flow pattern which is a solution to the original problem, we 
must subtract all the artificial flows of 5 along Aja, Aza, B,b; and Byb, 
from the flow pattern given above. This gives the following solution, in 
which we have indicated the direction of flow along each edge by an 
arrow. The original supplies and demands are shown next to the vertices. 


The total transportation cost is obtained by multiplying the value of the 
flow along each edge by the original cost associated with that edge. Thus 
the total transportation cost is 


(2x1)+(2x3)4+(1 x 4) + (2x 4) = 20. 


es 


Solutions to the problems 


Solution 1.1 
(a) The bipartite graph representing this situation is given below. 


applicants jobs 
1 A 
2 B 
3 C 
4 D 


(b) The only matching containing four edges is the following. 


applicants jobs 


1 A 
2 B 
3 fe 
4 D 


This is a maximum matching, since there are clearly no matchings 
with more than four edges. The matching shows how all four 
applicants can be assigned to jobs for which they are qualified. 


Solution 1.2 
(a) 
numberof Is marriage 
number women condition 

subset of men womenknown ~~ ofmen known satisfied? 
{A,, Az, A3} B,, Bo, B3, Bg 3 4 yes 
{A,, Ao} B,, Bz, Bg Z 3 yes 
{A,, A3} B,, Bo, B3, Bg 2 4 yes 
{A>, As} By, Bz, By 2 5 yes 
{A;} B,, B3, Bg 1 3 yes 
{A>} B, 1 1 yes 
{A3} Bo, B3 1 2 yes 


The marriage condition is satisfied for every subset, so, by the 
matriage theorem, all three men can marry women they know. 


(b) The number of non-empty subsets of a set of size n is 2” — 1. For this 
problem, n = 3, so we need to consider 2° —1=7 subsets — this agrees 
with the number of entries in the table. 


(c) The bipartite graph representing this situation is shown below. A 
maximum matching of three edges (corresponding to all three men 


being able to marry women they know) is indicated by thick lines. Other matchings with three edges are 
possible. 
men women 
B, 
Ay 
By 
A, 
Bs 
A; 
By 
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Solution 1.3 


(a) The bipartite graph representing this situation is: 
men women 
Ay By 
Ay By 
A; Bs 
Ay By 
(b) The only subset of men for which the marriage condition is not 
satisfied is the subset {A,, A,, A,} which contains three men who 
collectively know only two women, By and B3. 
(c) Since we know from part (b) that the marriage condition is not 
satisfied, the marriage theorem tells us that not all of the four men 
can marry women they know. It follows that, if we can find a 
matching containing three edges, then this is a maximum matching. 
There are several such matchings, two of which are: 
men women men women 
A a A B, 
Ay By A, By 
x B3 Ag B3 
Ag By Ag B4 
Since a maximum matching contains three edges, the maximum 
number of men who can marry women they know is three. 
Solution 1.4 
(a) The marriage condition is clearly not satisfied, since the five men 
between them know only four women. 
(b) The bipartite graph representing this situation is: 
men women 
Ay B 1 
As es 
A; Bs 
Ag By 
As 
(c) The subsets of men for which the marriage condition is not satisfied 
are listed below, together with the value of the expression 
p+(n —m) in the modified marriage theorem. 
number of 
number of women 
subset of men womenknown menm knownp p+(n-m) 
{Az, Az, Ag} By, Bs 3 2 3 
{A;, Az, A3, Ag} B,, Bo, B; 4 Ks 4 
{A,, Az, A3,A4,A5} By, Bo, B3, By 5 4 4 
(d) From the above table, we see that the minimum value of the 


expression p + (n — m) in the modified marriage theorem is 4. By the 


theorem, the largest number of men who can marry women they know 
is 4. 


All the subsets of men for which p + (n — m) is equal to its minimum 
value of 4 contain the men A», A3 and A,. It follows that if we add 
one additional woman Bs who knows any one of these men, then the 
marriage condition is satisfied for all possible subsets. Therefore, if 
we remove any one of A,, A3 and Ay, then the remaining four men will 
be able to marry women they know from the original four women. The 
list of subsets of four men who can all marry women they know is 
therefore: 


{A,,A3, Aq, As}, {A1,A2,Ag, Ast, {A1, A2,A3, As}. 


Solution 2.1 
(a) The path x3 y2 x2 y; is not an alternating path, since the initial vertex 
x, and the final vertex y,; are both incident with edges in M. 


(b) The path y3 x3 yz X2 y; X; is an alternating path. 


(c) The path x3 y4 X4 y3 is not an alternating path, since the initial vertex 
x3 is incident with an edge in M. 


Solution 2.2 
There are six alternating paths: 


X3 Y3 X2 Yj X Yr, X3 YW4 X4 Ys, X5 YW X4 Ys, X3 Yq X4 Y3 X2Y1 X%1 Yn, 
X5 Yq X4 Y3 X2 Y1 X1 YD, X3 Y2. 


Solution 2.3 


Let us start with the matching shown by thick lines in the bipartite graph 
below. 


Applying the labelling procedure, we obtain the following labelled 
graph. 


(Yo) xy Yi, (x) 
(y3) Xz Yq (Xp) 

x3 Y3 (x4) 
(+) x4 Y4 


Breakthrough occurs at y;. Tracing back through the labels from y;, we 
obtain the alternating path x4 y3 x2 y2 x;y. 


We thus obtain the improved matching shown on the left below and a 
corresponding maximum flow of value 4 shown in the network on the right. 
This is the flow which we would obtain by applying the maximum flow 
algorithm to the given network. 


»4 : 

. h ! 4,1 Val 

X2 Yo row 
Ky > 

X3 ¥3 ee a 

X4 Y4 Pili ia 


The number of edges in a maximum matching in the graph is equal to the 
capacity of a minimum cut in the network. 
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Solution 3.1 


The lowest cost in the cost matrix is 1, so we set p = 1 and subtract 1 from all 
the entries in the cost matrix. This gives the following cost matrix and 
partial graph, corresponding to p = 1. 


x4 Z 0 4 9 x4 OY, 
146s. Tas a 


X3 7 7 4 1 x3@ OY 


X4 a 1 x4@ OY, 


The partial graph has only one edge, so a maximum matching comprises 
just this edge, shown by a thick line. 


A complete assignment has not been found, so we repeat the process. The 
lowest non-zero cost is 1, so we subtract 1 from all the non-zero entries in 
the cost matrix. This gives the following cost matrix and partial graph, 
corresponding to p = 2. 


x4 Oy; 
oe 
lie 2 
a oe eee Fags 


We can now include the edges x2y3 and x3y, in the matching. 


A complete assignment has not been found, so we repeat the process, giving 
the following cost matrix and partial graph, corresponding to p = 3. 


XQ Y2 
/ \N 
X4 5 S y 0 x4 Y4 


There are no alternating paths in the partial graph, so the matching 
cannot be improved. 


A complete assignment is still not possible, so we repeat the process, 
giving the following cost matrix and partial graph, corresponding to 
p=4. 
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There are no alternating paths in the partial graph, so the matching 
cannot be improved. 


A complete assignment is still not possible, so we repeat the process once 
more, giving the following cost matrix and partial graph, corresponding to 
p=od. 


xy 0 G- @ 5 x4 Yj 
Xo \ a ee 0 Xp Y2 
X3 . ae 0 0) x3 ¥3 
x so ee 0 X4 Y4 


There are two alternating paths in the partial graph: 
XgY3XoYoXYr aNd X4Y4xX3¥3%QYoX1Y1 
giving the matchings 
X1 Yr, X22, X3Ya, XaYs aNd XYy,X2Y2, X3Y3, Xa Ya, 
respectively. These are both optimum assignments. 


A complete assignment has been achieved, and the minimum time for 
completion of an item is 5 minutes (the final value of p). 


Solution 4.1 


The application of the algorithm is outlined below. For convenience, we 
have omitted the circled numbers representing supplies and demands from 
the cost matrices. 


We assign weights to the vertices, and construct the first revised cost 
matrix and first partial graph. 


| ie | 0 


ba (2) 
Parts A and B 


Following the procedure of Steps 1-5, we find a maximum flow in the 
above partial graph. We thus obtain the following flow pattern. 


The edges carrying a flow are shown as thick lines. We have indicated 
the amounts of flow, and have subtracted these amounts from the supplies 
and demands at the appropriate vertices. 
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We return to Part A. 
Part A 


We carry out the labelling procedure. We label vertices A; and A3 with 
(+), and then label vertex b, with (A;) and vertex b3 with (A3). 


(*) (2) Ay 


No more labelling is possible, and breakthrough has not been achieved, so 
we go to Part C. 


Part C 


We find that 6 = min (1, 2, 3) = 1. 
We obtain the following revised cost matrix and revised partial graph. 
0 -1 0 0 


Part A 


We carry out the labelling procedure and obtain the following labels on 
the partial graph. 


Breakthrough is achieved at b,, and also at by, so we go to Part B. 
Part B 


The procedure of Steps 6 and 7 results in a flow of 1 being sent from A, to b1, 
and a flow of 1 being sent from A; to by. This results in the following flow 
pattern. 
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Part A 


We carry out the labelling procedure, and thus obtain the following labels 
on the partial graph. 


Breakthrough is achieved at b,, so we go to Part B. 
Part B 


Tracing back through the labels from b4, we obtain the flow-augmenting 
path Ab, Apb,. 


The even edge A,b, has a flow of 4. Thus the maximum possible backward 
flow is 4. 


The available supply at A, is 1. 

The required demand at by, is 1. 

The maximum flow F is therefore min (4, 1, 1) = 1. 

We send a flow of 1 along each of the odd edges A, b, and Ay by. 
We decrease the flow on the even edge Ap b; from 4 to 3. 


We thus obtain the following flow pattern. The original supplies and 
demands are shown next to the vertices. 


This flow pattern is the required minimum-cost solution. 


The total transportation cost is obtained by multiplying the value of the 
flow along each edge by the original cost associated with that edge. Thus 
the total transportation cost is 


(2x5)+(3x4)+(3x2)+(1x2)+(5x4) + (1x 4) = 54. 
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